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Topological superconductors are prime candidates for the implementation of topological-quantum- 
computation ideas because they can support non-Abelian excitations like Majorana fermions. We go 
beyond the low-energy effective- model descriptions of Majorana bound states (MBSs), to derive non- 
equilibrium transport properties of wire geometries of these systems in the presence of arbitrarily 
large applied voltages. Our approach involves quantum Langevin equations and non-equilibrium 
Green's functions. By virtue of a full microscopic calculation we are able to model the tunnel 
coupling between the superconducting wire and the metallic leads realistically; study the role of 
high-energy non-topological excitations; predict how the behavior compares for increasing number 
of odd vs. even number of sites; and study the evolution across the topological quantum phase 
transition (QPT). We find that the normalized spectral weight in the MBSs can be remarkably 
large and goes to zero continuously at the topological QPT. Our results have concrete implications 
for the experimental search and study of MBSs. 



Topological quantum computing [1 is interesting not 
only because of the technological goal but also because it 
requires the realization and manipulation of new types of 
quasiparticle excitations that have never been observed 
before and the question of their actual existence tests 
the limits of our understanding of the quantum world at 
a fundamental level. The simplest proposals along this 
line involve the realization of Majorana fermions (MFs). 

One of the pressing theoretical problems in this con- 
text is to sort out the experimental signatures of MFs. 
Several proposals for their realization in (effective) one- 
dimensional (ID) topological (p-wave) superconductors 
have been put forward recently [2J. A natural way to 
probe superconducting (SC) nanowires is via transport 
measurements and STM spectroscopy. Such transport 
studies of mid-gap features in the spectrum of unconven- 
tional superconductors have consistently attracted con- 
siderable attention [3] and particularly so in the quasi- 
one-dimensional case of (possibly p-wave) organic super- 
conductors [H [5]. The renewed interest in MBSs has 
given rise to more detailed studies including some excit- 
ing recent experimental claims of MBS detection [6|. 

Important insights have been gained from past theo- 
retical works [7HH] that study low-energy effective models 
consisting of the MBSs and an effective coupling to and 
between them as proposed in Ref. |7] These works have 
made specific predictions for detecting MBSs in exper- 
iments. For instance, non-equilibrium electrical trans- 
port between a metallic probe and a collection of coupled 
MBSs has been derived [7H9] and the differential electri- 
cal conductance calculated in the limit where the applied 
voltage (eV) is much smaller than the SC energy gap. It 
is crucial to go beyond low-energy approximations and 
calculate the full current-voltage characteristics includ- 
ing the contribution of the non-topological higher-energy 
excitations. It is also important to check the regime of 
validity of these low-energy models and even better have 
a full microscopic calculation to predict and interpret ex- 



perimental results in realistic systems. 

While different proposals for realizing MBSs in ID dif- 
fer in details 0, all essential features of realistic sys- 
tems are expected to be minimally described by the so 
called Kitaev chain (linked to the early theoretical pro- 
posals for topological quantum computation) [10 . Wc 
thus study non-equilibrium electrical transport in Ki- 
taev's tight-binding model of a single-channel ID quan- 
tum wire with triplet (p-wave) SC pairing between elec- 
trons with the same spin orientation (equal-spin pair- 
ing). The model's spectrum includes two MBSs at the 
two ends of the wire, which are expected to exhibit non- 
Abelian braiding statistics (in wire networks [H]), are 
separated from the non-topological complex excitations 
by an energy gap and have been proposed as a basis for a 
decoherence-free topological quantum memory. In view 
of theoretically capturing the physics of realistic exper- 
iments, we are interested in studying non-equilibrium 
transport through the ID quantum wire in an open- 
quantum-system geometry, i.e., the wire is connected to 
external baths which incorporate both decoherence and 
dissipation. This is a minimal setup that can capture the 
essential features of recent proposals for realizing p-wave 
SC wires |2J. We develop a novel extension of the quan- 
tum Langevin Equations and Green's Function (LEGF) 
approach [12] to SC wires with BCS-type pairing. LEGF 
provides us with an appropriate exact method to de- 
termine the electrical current in the SC wire with arbi- 
trary values of eV as compared with the strengths of the 
contacts and the SC pairing amplitude. Although this 
method has proven quite useful to investigate electrical 
[I2TH4] as well as thermal [15] non-equilibrium transport 
in non-interacting, ordered and disordered quantum sys- 
tems connected to baths, in this letter we will focus only 
on the electrical transport. 

The Hamiltonian of the full system is given by H = 
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The term Hw denotes the Hamiltonian of the ID wire 
and "Hg that of the p th bath, while the coupling between 
them is W^, B . Here a/ and dP a denote standard fermion 
operators on the SC wire and on the p th bath, respec- 
tively. The p-wave BCS-pairing amplitude of the wire is 
A and the parameter v can be tuned to drive the SC wire 
in or out of a topological phase (cf. Refs. [10] and [16]) . 
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FIG. 1. Schematic depiction of a Kitaev chain in terms of MFs 
(filled circles) for v = 0. Two generic cases with odd (N = 3) 
and even (N = 4) number of sites are shown side by side. The 
solid and dashed lines correspond to |A| ±7 matrix elements, 
respectively. The two dotted lines at each end represent equal- 
strength contacts with the respective lead. Since for v = 
there is no matrix element between the two MFs within each 
site, one can see the existence of two independent Majorana 
(conduction) channels with N MFs each. Notice those two 
channels are always mirror-equivalent for odd N and contain 
an exact zero mode each. For |A| = 7 and any odd or even 
N, the dashed lines are absent and there are two zero-energy 
MBSs localized one at each end and decoupled from the rest 
of the chain. 

Next we introduce a local re-parametrization for 
the electron operators on the SC wire in terms of 
MFs lift, ai = e- ie l 2 {c 2 i- x + ic 2l )/2, with c m = 
c m and {c m ,c m '} = 2d m>m ' . The wire Hamilto- 
nian in the MF basis (see Fig. [I]) is given by = 

(i/2) [(-!>) E^iQw-K* + E ; =i {(|A| + 7) c 2l c 2l+1 + 
(|A| — 7)c 2 /_iC 2 ;+2}] ■ The two MFs at any one end of 
the wire are connected to the end site of the bath with 
equal-strength coupling. (Note that the contact coupling 
to the inner MF of the end site of the wire is outside 
the scope of all previous works, as they deal either with 
effective MBS models [THS] or a continuum representa- 
tion [3].) We assume that for times t < to, the wire is 
disconnected from the baths (each in equilibrium at a 
specified temperature T p and chemical potential fi p ). At 
time to we connect both baths to the ID SC wire and we 



are interested in the steady-state properties, which can 
be derived by integrating out the bath operators using 
the quantum LEGF formalism and considering the limit 
t — > —00. We get the following solution 
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where c m (oj) = (l/2n) dt e lLUt c m (t) and same for 
fj p (uj). The single-particle retarded matrix Green's func- 
tion of the p th bath is given by g p+ {t) = -i6>(t)e- m B*/ n 
and gfjt (w) = dt e^g^t (t) . Its first diagonal ele- 
ment (in the MF 'site' basis) is used to define E+(cj) — 
(7p/^) 2 3ii~( w )j which corresponds to dissipation due to 
coupling with the bath. The 'noise' in the leads, 
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depends on the bath's initial distribution and obeys |13] 

= r pM fi^' Vp> T p) <K W ^ w ') <W( 4 ) 

where f{cj,fj,,T) is the Fermi distribution function. The 
above identity is a form of fluctuation-dissipation rela- 
tion. We use the definition T p = — Im[E+]/7r = {l' p /K) 2 p p 
where p p {ui) is the local density of states at the first site 
on the p th bath. Finally, G + (w) is the Green's function 
of the full system (wire plus baths) [15] . 

The outward current between the wire and the p th 
baths is given by 
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where (...) denotes averaging over noise like in Eq. Q. 
For the sake of concreteness and simplicity of presen- 
tation, we will consider the 'symmetric' case for which 
the two baths to are similar (7^ = 77?) and connected 
with the SC wire by identical contacts (7^ = 7^ or 
T L = T R = r) while fi L = fl/2 = —fx R and T L = T = T R . 
In this case the steady-state currents are equal up 
to a sign and 
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We have checked that the expression in Eq. (|6| reduces 
to the known result for the case of a normal metallic 
wire (A = 0) where T{ui) is the transmission coefficient 
of a single electron through the wire. For a SC wire, 
T(u>) that contains contributions to transport coming 
from both electrons and holes, can be simply interpreted 
in the symmetric case via the zero-temperature differen- 
tial conductance (with A = sV, e = 1) 
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Moreover, for v = 0, T(ui) is an even function and thus 
proportional to dl/dV. 

At v = and 7 = A, there exist two exactly zero- 
energy localized MBSs (decoupled from rest of the wire) 
and the corresponding zero-temperature dl/dV has a 
peak at zero bias and two more peaks in dl/dV at 
eV = ±(7 + A) coming from all the degenerate high- 
energy excitations at these parameter sets. As can be 
seen from Fig. [2ja-f), the dl/dV shows distinct behav- 
iors for finite TV depending on whether TV is odd (panels 
e,f) or even (panels a-d) once 7 ^ A. We have developed 
a physical understanding of different behaviors by closely 
looking at the behavior of a finite Kitaev chain with odd 
vs. even TV (see Fig. [T]). For odd TV, there are TV pairs of 
degenerate modes including a pair of zero-energy modes 
thereby showing TV peaks (including one at V = 0). For 
even TV, all modes, including the ones at zero energy are 
split and the zero-temperature dl/dV has a dip at zero- 
bias. However, the splitting in the MBSs decreases with 
increasing wire length, and finally goes to zero in the 
long-wire limit as exp(— aN) with a ~ (A + 7)/|A — 7I 
for A ~ 7. (This is in contrast to the result of an ef- 
fective model of a chain of coupled MBSs, for which the 
zero-bias conductance is zero for the case of two MBSs 
[jj].) Note that the splitting of the non-topological high- 
energy excitations survives for a long even-site wire as 
can be seen by comparing heights in Fig. [2ja-d) . 

For a finite v 1 the high-energy excitations are split for 
any even or odd number of sites |17j . Gating the wire 
(0 < \u\ < 2|7|) also splits any degenerate zero modes 
at V — for the finite TV but, remarkably, in the large 
TV limit, they regain their integrity and yield a zero-bias 
peak. (The splitting for finite- TV case can be understood 
by adding a matrix element between the two MFs within 
each site in Fig. [I]) Next, in the long- wire case, we closely 
study how dl/dV evolves as v is varied across the topo- 
logical QPT [i.e. the quantum phase transition from the 
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FIG. 2. dl/dV (in units of e 2 /h) vs. eV for 7p = 1, j' p = 
7 = 0.1 (p = L, R) and |A| = 0.3. Everywhere v — except 
for v = 0.1,0.25 in (g,h) respectively. Wire lengths are as 
follows: (a) N = 2, (b) TV = 6, (c) N = 16, (d) N = 100, and 
(e) N = 3, (f-h) TV = 101. 



topological < 2|7|) to non-topological > 2\j\) 
SC phase]. We observe that as \v\ approaches 2|7|, the 
gap shrinks and the zero-bias peak corresponding to the 
MBSs loses weight; at \v\ — 2\^\ the gap closes com- 
pletely; and for \v\ > 2\j\ the gap reopens with no zero- 
bias anomaly [see Fig. [2](h)] . 

In Fig. [3^a) we compare the I-V characteristics of the 
wire for different parameters (varying A). In the topolog- 
ical phase, I displays a sharp rise for small voltages due to 
the presence of the low-energy MBSs and then saturates 
with increasing voltage reflecting the gap in the energy- 
spectrum. The gap and hence the plateau is most sharply 
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FIG. 3. I-V characteristics for a symmetrically biased wire. 
N = 15, 7 P = 10% = 1 (p = L,R) and 7 = 0.3. In (a) A 
changes as 0, 0.1, 0.6 while v — 0; and in (b) v = 0.65. 



defined for |A| > |7|. The current increases again for 
even higher voltages, when non-topological excitations 
in the SC wire start to contribute to the transport, and 
finally saturates again when the applied voltage is larger 
than the band-width of the wire. For comparison, we 
plot in Fig. |3^b) the I-V of the topologically trivial SC 
phase along with the corresponding metallic phase with 
the same v. 

We find that the normalized spectral weight of the 
MBSs initially decreases with increasing wire length [see 
Fig. Qa)] but quickly saturates to a length-independent 
finite value (cf. with the ideas of electronic teleportation 
in these systems [18]) which depends on 7 and A but is 
independent of the strength of the contacts. As can be 
seen in Fig. Hb), the maximum spectral weight for long 
wires is 50% (a remarkably large fraction carried by a sin- 
gle pair of states) and is achieved when 7 = A and v = 
when the two MBSs have exactly zero energy, are com- 
pletely localized (one at each end) and decoupled from 
the rest of the wire (see Fig. [I]). The contribution of the 
MBSs (which can be interpreted as the order parameter 
of the topological phase) decreases for increasing v and 
at the topological QPT goes to zero in the way shown in 

Fig.gb). 
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FIG. 4. Normalized spectral weight of the MBSs for 7 P = 
10 7; = 1 (p= L,R), and (a) A = 0.3 or (b) 7 = 0.3. 

If we depart from the symmetric case considered above 
by setting fiR = 0, we find that j'l gives an I-V similar 
to that of the symmetric case while j'r becomes zero (in 
agreement with Ref . [!§]) . We have also verified that mod- 



erate asymmetry between the left and right baths or the 
respective contacts does not significantly affect the re- 
sults for the current as compared to the symmetric case 
and is thus meaningful to relate these to existing or fu- 
ture measurements. 

In conclusion, by carrying out an exact microscopic 
calculation of non-equilibrium electrical transport in a 
p-wave (Kitaev) chain connected to baths, we have been 
able to go beyond the previous works using either effec- 
tive models of MBSs [7HS] or a continuum representation 
such as the Bogoliubov-deGennes formalism |4j. First, we 
have captured the interplay and differences between Ma- 
jorana and regular complex fermions (like the Bogoliubov 
excitations in the superconductor) and overcome the lim- 
itation of having voltages smaller than the SC gap. Sec- 
ond, we have modeled the contacts more realistically and 
with greater detail (going beyond the scope of scattering- 
formalism approaches; cf. Refs.0]and[Hl). Third, we have 
shown how the microscopic properties and correspond- 
ingly the differential conductance varies with the length 
of the SC quantum wire for different parameter regimes. 
Fourth, we have compared and contrasted the behaviors 
when the number of sites in the wire is even vs. odd. 
Fifth, we have calculated the normalized spectral-weight 
contribution of the MBSs to the total differential conduc- 
tance and thereby provided a starting point to meaning- 
fully compare the contributions to the transport coming 
from MBSs vs. other excitations in realistic experiments. 
And finally, we have predicted how the differential con- 
ductance changes across the topological QPT and shown 
the way in which the spectral-weight of the MBSs con- 
tinuously goes to zero at the QPT. 

As shown in this letter, all these points not only make 
concrete experimental predictions, but also have strong 
implications for the scope and validity of the low-energy 
effective models. Key features of the experimental dl/dV 
curves E] are in agreement with our results, in partic- 
ular the overall shape of the zero-bias conductance peak 
and the pronounced dips on its sides (which reflect the 
plateaux shown in Fig. [3]). We suggest systematic mea- 
surements of the differential conductance based on our 
predictions and, in particular, the analysis of the spec- 
tral weight of the zero-bias conductance peak (such study 
would require a careful material-dependent protocol to 
take into account the effect of higher bands and dimen- 
sionality when defining the fractional spectral weight) as 
a way to strengthen and advance the experimental efforts 
to observe MBSs. 
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